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Cross-correlation of WMAP7 and the WISE Full Data Release 
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ABSTRACT 

We measured the cross-correlation of the Wilkinson Microwave Anisotropy Probe (WMAP) 7 
year temperature map and the full sky data release of the Wide-field Infrared Survey Explorer 
(WISE) galaxy map. Using careful mapmaking and masking techniques we find a positive 
cross-correlation signal. The results are fully consistent with a ACDM Universe, although not 
statistically significant. Our findings are robust against changing the galactic latitude cut from 
\b\ > 10 to \b\ > 20 and no color dependence was detected when we used WMAP Q, V or W 
maps. We confirm higher significance correlations found in the preliminary data release. The 
change in significance is consistent with cosmic variance. 
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1 INTRODUCTION 

Cosmological supernovae measurements and Cosmic Microwave 
Background (CMB) fluctuations support cosmological models in 
which the cosmic energy density is dominated by Dark Energ y 
(DE) at the present epoch djarosik et all 1201 it iRiess et al.lll998T) . 
In such theories the current accelerating expansion and the decay 
of gravitational potentials are predicted. Therefore, the presence of 
DE is manifested in both geometrical and dynamical forms. 

Dark Energy comes to dominate the energy density at late 
times, z < 2, and so the primordial fluctuations in the CMB alone 
do not provide a sensitive probe. However, DE may leave a signal 
in the secondary anisotropies that are imprinted on the microwave 
background radiation. The Integrated Sachs -Wolfe effect (ISW) 
dSachs & Woife1ll967l : iRees & Sciamalll968l) is an example of a 
secondary anisotropy: CMB photons passing through a changing 
graviational potential become slightly hotter or colder. In a flat and 
matter-dominated Universe the potential is constant on large scales 
thus gravitational blueshifts and redshifts cancel along the photon 
path. However, in a Universe dominated by DE there is a net energy 
difference between entering and leaving a potential well due to the 
decay. Thus, the detection of the linear ISW effect provides direct 
evidence for dark energy in the ACDM model. Furthermore, alter- 
native gravity models provide predictions for the ISW effect and 
may b e directly tested with ISW observations dGiannantonio et all 
l20ich . 

The ISW signal may be detected through cross-correlation of 
large-scale structure surveys with the CMB temperature maps. The 
correlation is weak, generally less than a 1 fiK signal is expected, 
orders of magnitude below the primary fluctuations. Futhermore, 
the ISW effect is strongest on large angular scales where cosmic 



variance is also large, making the measurement even more cumber- 
some. 

Several measurements have been performed to uncover 
the ISW signal: positive cross-correlations were measured us- 
ing galaxy dat a from the Sloan Digital Sky Survey ( SPSS) 
and WMAP dFosalba et alj 120031: JPadmanabhan et alj 120051 ; 
lGranettetal.ll2008ll2009l:|PaDaietal.ll201ll) . Other successful at- 
tem pts were |Fosa lba & Gaztana gal (120041) base d on APM galax- 

using radio data 



tem pts were rosalba & 
ies. Wolta et alj <2004) 



Raccanelli ct 



aL N2008h 



and iBoughn & Crittenden! d2004d lbl) in which the har d X-ra y 
background was invest igated. Besides, Afshordi et al. I d2004h : 
iRassat et alj d2007l) and iFrancis & PeacockT i 2010h used infrared 



galaxy samples to characterize the ISW signal. The typical ISW 
significance in the papers above is around 2-3cr. Comp rehensive 
studie s usi ng combinations of data sets were carried out bv lHo et all 
d2008l) and lGiannantonio et alj d2008l.l2012h . 

The Wide-field Infrared Survey Explorer (WISE) all-sky sur- 
vey is an attractive dataset for ISW studies. The survey effectively 
probes low redshift z < 0.3 with a high source density. Us- 
ing the prelimina ry data release (PDR) covering 10000 sqr deg. 
lGotoetal.ld2012l) cross-correlated a WISE galaxy sample with the 
Cosmic Microwave Background, finding a 3<r detection, although 
with three times the amplitude expected in ACDM. In this paper, 
we rexamine this finding using the full-sky data release (FDR) of 
the WISE survey and the WMAP 7-year dataset. 

The structure of this paper is as follows. In Section 2 we de- 
scribe the data we used in particular. Section 3 describes our meth- 
ods including the theoretical expectations, simulations and mea- 
surements. Finally, in Section 4 the statistical significances are pre- 
sented and systematic effects are discussed. 
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Figure 1. WISE allsky galaxy sample, together with our mask including the 
stripes and WMAP's mask area. 



2 CMB DATA AND GALAXY MAP 

We used the best achieveable versions of the CMB data products 
and focused on the reliability of our new galaxy sample. In this 
section we describe our map-making and masking procedures. 

2.1 CMB map 

The 7-year WMAP t emperature data we re downloaded from the 
LAMBDA websiteQ Jjarosik et af1l201 ll) . WMAP data is affected 
by noise and contaminations both from point sources and the Milky 
Way. Among all, the Q, V and W maps have the least galactic con- 
tamination. We used the foreground reduced version of these maps 
and the CMB Extended Temperature Mask was chosen to avoid 

contaminations. 

Using HEALPIX l lGorski et al.l |2005|) . NSIDE=128 repix- 
elized versions of the maps were created. Galactic foregrounds and 
known point sources are quite surely excluded and 71% of the sky 
is unmasked. 

2.2 WISE galaxy map 

The density map of the gal axies was prepared using the full data re- 
lease of the WISE project jWright et al.ll201(]|) . The WISE satellite 
surveyed the sky at four different wavelengths: 3.4, 4.6, 12 and 22 
/im. We used different ban ds to separate stars from galaxies using 
color-color plots. Following lGoto et al.ld2012h we select sources to 
a flux limit of Wl ^ 15.2 mag to hav e a uniform dataset. 

According to iGoto et alj J2012l) the majority of stars near to 
the galactic plane have a W3.4 — W4.6 ^ 0.2 mag color. Moreover, 
it was found that a Wa.6 ~ W12 ^ 2.9 mag selection reduced the 
stellar contamination. We confirm these findings in the FDR and 
followed the same procedure for star-galaxy separation. 

Our galaxy sample exhibit s stripe-like overd ensities on the 
map in several locations. While IGoto et all d2012h applied hand- 
made cutouts in their mask to exclude regions with unusually high 
number counts, we understand that the stripe-like features origi- 
nated from the observational strategy of WISE and the position of 
the moon. We realized that the moon-contamination flag may be 
used to properly mask these regions. We forced out pixels in which 
the 'moonlev' flag is higher than 3 in at least one of the bands. This 

1 http://lambda.gsfc.nasa.gov/ 
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Figure 2. Gradients in galaxy density can be detected as a function of galac- 
tic latitude b. The effect is less significant if objects with warning flags are 
included in the sample, since the density drop is 4% at |6| = 20 com- 
pared to its value at |b| = 50. 



Table 1. Mask properties of the WISE galaxy sample. The first 3 rows cor- 
respond to allsky data with b galactic latitude cuts. The last mask was made 
using only the preliminary survey area and an additional |6| > 10 cut. 



Mask 


Areafdeg 2 ] 


jynoflag 
gal 


jystarcorr 
gal 


N gal 


|6| > 10 
6 > 15 
6 > 20 


26443 
25248 
23167 


1662995 
1617745 
1521755 


2016563 
1948454 
1814156 


2234370 
2129013 
1945517 


161 > 10 


10967 


782502 


976209 


1057073 



means that the fraction of the used image frames suffered by moon- 
contamination is higher than 30%. Masking regions based upon the 
moon contamination flag effectively removes the stripe pattern. We 
cannot address any further residual effects of the moon contamina- 
tion outside our mask area. We have found that with a more con- 
servative magnitude limit of Wl < 14.9 overdensities along stripes 
are reduced in width but do not disappear. 

Sources in WISE may also be contaminated by artefacts in- 
cluding the halos of bright stars, ghost images and diffraction 
spikes and special data-quality flags exist to handle these problems. 
If the dataset is filtered using these flags we possibly lose real galax- 
ies. However, if we do not use flags to create a conservative catalog 
then stars or galaxies w ith insufficient parameters can appear in the 
data. IGoto et alj d2012h used additional pixels in their mask to ex- 
clude regions where the abundance of these unreliable objects is 
high, but did not filter the whole galaxy sample. We investigated 
both cases and we have found this choice is important especially 
on large scales. 

We find a gradient in the galaxy density with Galactic latitude 
with fewer galaxies ne ar the Galactic pl ane. An empirical correc- 
tion was developed bv lGoto et alj | |2012|) in which a mean correc- 
tion is computed in galactic latitude bins to artificially 'flatten' the 
distribution at |fc| < 20. 

We attribute the gradient to stars near the Galactic plane mask- 
ing background galaxies. The problem is made more severe due 
to the broad point spread function of WISE (6-12 arcsec). We use 
the Tycho2 star catalogue which reaches a depth of V < 13 mag 
dH0g et alj|2000h to measure the survey area lost around each star. 
We calibrated a mean relation between V magnitude and star halo 
radius R for WISE R = 9.52 - 0.741/ for R in arcminutes. Any 
detected sources within this radius of a Tycho star is removed. We 
then construct a map of the lost area by summing the area attributed 
to stars in each Healpix pixel. This map is then used to normalise 
the galaxy counts. Figure|2]shows the result of this correction. 
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Apparently, the density gradient does not disappear, so only 
a modest correction is possible with our method. Interestingly, we 
find that the gradient is higher when flagged sources are not in- 
cluded in the sample. The exact source of the gradient is unex- 
plored, but we show that this effect did not mess up our measure- 
ments, our findings are robust. 

As described, regions nearby the plane of the Milky Way are 
potentially contaminated, and we consider the most appropriate so- 
lution to use \b\ > 20 regions. To perform tests with only the pre- 
liminary sky coverage area a mask was created using the area cov- 
ered by the preliminary survey. 

2.3 Redshift distribution 

In order to calculate a theoretical ISW expectation, redshift in- 
formation is needed. Since WISE is a photometric survey with- 
out spectroscopy, the selected galaxies were cross-id entified with 
source s from GAMA (Galaxy and Mass Assembly, iDriver et al.l 
(2011)) sample that has spectroscopic redshift of ~ 200000 galax- 
ies. Using the overlapping part of the two surveys we have found 
a pair for 82% of the galaxies with a 3" matching radius. We esti- 
mated an accidental matching rate of 0.1% for this analysis using 
random points with WISE density. The matched sample has a me- 
dian redshift z ~ 0.15. The obtained approximate redshift distribu- 
tion provides a basis to calculate a theoretical cross-power spectrum 
in this redshift range. 



3 RESULTS 

In this section we discuss the results using the CMB and galaxy 
datasets. We also elaborate on the most important theoretical and 
simulated considerations related to ISW detection and further anal- 
ysis. 

3.1 WISE-WMAP cross-correlation 

We calculate the cross power spectrum using a fast quadratic 
esti mator SpICE ( Spatial l y Inhomogen ous Correlation Estima- 
tor, ISzapudi et af] J200ll 120051, 1200 lh ). The individual band 
powers are binned logarithmically, the boundaries are I — 
6, 8, 11, 16, 22, 31, 44, 61 and 87 this means the first band stands 
for 1 — 6,7 etc. With this choice we avoid the lowest / range, where 
c osmic variance is m eaningful and it is easier to compare the results 
to lGotoetal1 ( l2012l) . who used the same bins. Our measurement is 
introduced on Figure [3] 
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Figure 3. Cross-correlation power spectra of WMAP7-WISE datasets us- 
ing our |6| > 20 mask, together with theoretical expectations for ACDM 
cosmology and zero ISW detection. Magenta errorbars were computed us- 
ing 1000 simulations and the diagonal elements of the covariance matrix. 
The light gray errorbars represent the variance of the Ci FDR — Ci PDR 
differences in each bin to check if the change is consistent with cosmic 
variance. 



function. Independent determination of b g is not possible in lin- 
ear theory, because erg acts to renormalize the power spectrum, 
Ci oc ( ba$ ) 2 . Thus we fit only for b g and keep as fixed. CosmoPy 
and CAMB rl were used to gen erate nonlinear matter power spectra 
with Halofit Tsmith et alj2003h at the median redshift of our galaxy 
sample. 

We measure the galaxy-galaxy power spectrum with SpICE. 
The measurement is affected by Poissonian shot noise that has a 
form of 1/JV, where N is the mean galaxy count per steradian. 
Actually, the impact is negligible, less than 10% at the maximum 
I we used and less significant at larger scales where we expect to 
measure ISW. However, to be precise we substracted the noise and 
the amplitude of the theoretical model curve was fitted on the 6 < 
I < 100 interval, i.e. angular scales down to ~2°. Our result is 
b g = ba 8 = 1.04 ± 0.05. 

To bring the b g parameter into use consider now the expression 
of the theoretical ISW signal. The cross-spectrum of a galaxy map 
and the CMB is given by 

n 3T l 6 • TcMB^lmHo 2 f , 2 D 

Cf =b g / dkk P k - 

J drj,(kr)k~ 2 d ( 1 +2 Dl{z) J dr'Mkr^^r'y 2 (2) 

where D 1 (z) is the linear growth facto r, the numerica l result of this 
expression depends on the cosmology dCooravll2002l) . 



3.2 Theory 

We derive the expected correlations and galax y bias using WMAP 7 
best-fit A CDM cosmologic a l para meters djarosik et al.l 1201 ll) . 
Following iFrancis & Peacock! i hoiOl) a linear bias relation is con- 
sidered to couple galaxy and matter overdenisities, S g — b5 m - The 
two-dimensional projection of the 3D galaxy auto-correlation is 
given by 



C gg (l) 



dk k 2 P(k) 



drr <j)(r) ji(kr) 



(1) 



where 4>{r) oc —^p-^fr ' s a com oving coordinate with a nor- 
malization relation J 4>{r)r 2 dr = 1, and ji is a spherical Bessel 



3.3 Simulations 

We simulated 1000 random CMB skies with AC DAI cosmologi- 
cal parameters using Healpix synf ast to cross-correlate with our 
WISE galaxy density map. The power spectrum was calculated and 
binned into the 8 spectral bins given above. The covariance matrix 
estimated from these measurements is shown in Figure [4] Neigh- 
boring bins are anti-correlated typically by 10%. The diagonal ele- 
ments were used to calculate the errorsbars shown on Figure [3] 



2 http : / /www .ifa.hawaii. edu/cosmopy/ 

3 http://camb.info/ 
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Table 2. Significance properties of our results are shown, x" 
table are less than theoretically expected, this indicates that the error bars 
are possibly overestimated. The same was reported by Francis & Peacock 
and Rassat et al. in their work. We performed Monte Carlo runs with 1000 
and 3000 trials, but the covariance was robust. More over, the analytic esti- 
mate of the covariances and errors (jcabre et al. 20071) agrees to 5% with the 
MC outcomes. 



Mask 


ISW Model 


x 2 


Ax 2 


Amplitude 


S/N 




Zero 


3.07 








|6| > 10 


Best-fit 


2.20 


0.87 


0.8 ±0.9 


0.9 




ACDM 


2.26 


0.81 








Zero 


2.71 








|6| > 15 


Best-fit 


2.13 


0.58 


0.7 ±0.9 


0.9 




ACDM 


2.27 


0.44 








Zero 


2.32 








|6| > 20 


Best-fit 


1.63 


0.69 


0.8 ±0.8 


1.0 




ACDM 


1.74 


0.58 






|6| > 10 


Zero 


5.64 








preliminary 


Best-fit 


2.91 


2.73 


2.3 ± 1.2 


1.9 


area only 


ACDM 


3.74 


1.90 







Figure 4. Covariance matrix of the logarithmic spectral bins from 1000 
random CMB simulations. Neighboring bands are not correlated in a good 
approximation. 



4 SIGNIFICANCE TESTS 

Again we follow iFrancis & Peacock! feoich . now to determine the 
significance of our ISW detection. Consistency of our results was 
investigated with three hypotheses: null-detection of ISW, regular 
ACDM model predictions and finally with a best-fit theoretical 
curve. Our statistics is based on the amplitude fit, we set the ampli- 
tude of the ACDM theoretical curve to 1.0 and to 0.0 in the zero 
ISW case. 



4.1 Statistical tools 

We evaluate a \ 2 statistic for each hypothesis which is the follow- 
ing: 

X 2 = ^rf l ^ 1 d J (3) 

ij 

where C is the covariance matrix and di = {Cf ata -Cf yvo ).Cf yvo 
can be given by Equation [2] assuming various models or it is zero 
in the null-ISW case. Index i labels the bins we use in the cross- 
spectrum. Moving forward one step we define the likelihood of a 
hypothesis below: 

Co,\C\- N/2 e-^ Tc ' lA (4) 

where N is the number of data points and with d vector that was 
constructed using di and C matrix. 

In fact, our tool to describe the detailed statistical properties 
of our test is —2 In (C1/C2) = A% 2 where the ratio of the two 
likelihoods is taken in a case of two different hypotheses and A\ 2 
is calculated. In general, A\ 2 > 3 is a strong evidence for a sig- 
nificant difference. 



4.2 Systematic effects 

Although we have taken care in source selection and masking we 
must address any residual systematic effects that may exist. Galac- 
tic dust or emission detected by WMAP could contribute to a corre- 
lation with WISE due to dust attenuation or the gradient in source 
density measured with Galactic latitude. However, we performed 
all the significance tests using Q, V and W foreground reduced 
CMB maps and the maximum relative difference in a given spec- 
tral bin was 1.2%. This fact lead us to the conclusion that there 
is no significant color dependence and effects from Galactic dust 
or emission must be minor. We also checked the results using the 
WMAP team's Temperature Mask or Extended Temperature Mask 
but no meaningful difference was found. 

On the other hand, we investigated several systematic effects 
related to our galaxy sample. Firstly, we analyzed different galac- 
tic latitude cuts. Our finding is robust, even if not significant, the 
results are summarized in Table [2] 

Next we considered the differences in the detection signifi- 
cance due to initial mapmaking. Our tests were applied to a map 
without flagged objects and it was obtained that the amplitude 
varies between 0.5 and 1.0 regardless of the starmask correction 
technique. 

We repeated our analysis with Wl ^ 14.9 mag limit to test 
the effects of faint sources. With this sample originating from a 
slightly different redshift distribution we found an increment in the 
ISW signal, but the errorbars were also high so the significance 
remained ~1.0cr. 

We also extended our cross-correlation analysis to 2 ^ / sj 5 
multipoles and very weak positive cross-correlation was measured 
with extremely high errorbars. The results increased the signifi- 
cance only slightly, however, they were sensitive to the galactic cut. 

In the light of the robustness of our results against different 
galactic cuts we argue that the stellar contamination is low or at 
least uniform in our galaxy sample and it does not affect the mea- 
surements. The upper limit is 18% from the WISE-GAMA match- 
ing but a similar deficit of optic al pairs was also reported using 
SDSS and WISE jYan et alj|2012h . These facts all indicate that the 
different selection criteria of the infrared and the optical bands are 
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responsible for the missing counterparts, rather than confusion with 
stars. In summary, the presence of stars is probably less than the un- 
paired fraction of our galaxy sample. 



5 DISCUSSION AND CONCLUSIONS 

Our principal aim was to produce the final ISW measur ement us- 
ing WISE, we compare our results with lGoto et al.ld2012h that used 
the PDR. We repeated our measurements with the WISE prelimi- 
nary data and largely reprod uced the individual Ci gT powers were 
found bv lGotoet ail d2012n . although we measured lower signifi- 
cance, except using alternative binning. We do not expect perfect 
agreement, given that the analysis was performed from the ground 
up. With the same bins our best fit amplitude for the PDR was 
2.5 ± 1.2 i.e. a 2.1cr detection. The result is consistent with our 
1. 9a finding on the preliminary part of the sky but using the new 
data. Using the full sky we measure an ISW significance of ~1.0<j. 
The change is fully consistent with possible cosmic variance, as 
illustrated by the light gray errorbars in Figure |3] 

With our enhanced mapmaking and better view of the WISE 
data we suppressed artefacts. Our mask is entirely based on the 
properties of the WISE object flags and many systematics were re- 
vealed. However, the signal decreased despite the improvements in 
our analysis methods. 

While some recent studies especiallv lGoto et"ai1d2012|) raised 
the possibility that the ISW correlations might be higher than 
ACDAI predictions, we conclude that the signal we found is 
consi s tent with ACDM and pre vious measurements dRassat et alj 
120071 : iFrancis & Peacockl |2010|) . Our analysis highlighted that 
higher ISW amplitude measurements on certain parts of the sky 
can be due to cosmic variance. 
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